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We investigated the effects of nonequilibrium and collision 
terms on the deconfinement phase transition of an expanding 
quark system in Friedberg-Lee model in relaxation time ap- 
proximation. By calculating the effective quark potential, the 
critical temperature of the phase transition is dominated by 
the mean field, while the collisions among quarks and mesons 

change the time structure of the phase transition significantly. 
PACS : 05.60. + w, 52.60. + h, 12.38.Mh 

I. INTRODUCTION 

It is generally accepted that the most interesting quan- 
tum chromodynamics (QCD) phase transition in hot and 
dense nuclear matter is the deconfinement phase tran- 
sition between normal nuclear matter and quark-gluon 
plasma (QGP), where quarks and gluons are no longer 
confined. The theoretical and experimental investigation 
of QGP is one of the most challenging problems in high 
energy physics. It is widely believed that QGP phase can 
be formed in ultra- relativistic heavy-ion collisions. 

Due to the difficulty of non-perturbative treatment in 
QCD, various models have been considered in the study 
of the phase transition among which Friedberg-Lee model 
[1], also referred as nontopological soliton model, has 
been widely discussed in the past two decades, see for 
instance [2,3] and refs. therein. In this model, the 
non-perturbative dynamics responsible for confinement 
in QCD is simulated in terms of a non-linear coupling 
to a scalar field a. It shows an intuitive mechanism for 
the deconfinement phase transition. In vacuum state, the 
physical value of a is large and the quark mass is more 
than 1 GeV, so that the effective heavy quarks have to 
be confined in hadron bags [1]. With increasing temper- 
ature and/or density of the system, the physical value of 
a and in turn the effective quark mass drops down, the 
thermodynamic motion leads to a deconfinement of the 
effective light quarks. 

Wilets and his cooperators [3,4] did a great deal of 
work on the properties of Friedberg-Lee model, mostly 
in mean field approximation and in vacuum state. It is 
proved very successful in describing the static proper- 
ties of the nucleon. During the past years, Friedberg-Lee 
model was extended to finite temperature and density to 
study deconfinement phase transition [5-8]. Similar to 
most of these investigations in the frame of finite tem- 



perature field theory, the temperature and density effect 
on the phase transition is based on the assumption of 
a thermalized plasma phase. While one can use vari- 
ous parameters and take different treatments, the criti- 
cal temperature of the deconfinement phase transition in 
the model is limited in the region of 80 to 120 MeV at 
zero chemical potential, much lower than the prediction 
in Lattice QCD [9]. 

Because of the estimated very short lifetime of the 
heavy ion collision zone, the highly excited particle sys- 
tem may spend a considerable fraction of its life in a non- 
thermalized, preequilibrium state. The dynamical tool 
to treat dissipative processes in heavy ion collisions and 
the approach to local thermal equilibrium is in principle 
nonequilibrium quantum transport theory. A relativistic 
and gauge covariant kinetic theory for quarks and par- 
tons has been derived [10], both in a classical framework 
[11,12] and as a quantum kinetic theory [13,14] based 
on Winger operators defined in eight-dimensional phase 
space [15]. To solve the quantum kinetic equations as 
initial problems, the transport and off-shell constraint 
hierarchies have been established [16-19] in the frame of 
equal-time Wigner operators. The properties of nucleon- 
nucleon collisions have been explored [20,21] in transport 
approach of Friedberg-Lee model. 

In the present paper, we consider an expanding non- 
equilibrium system with collision terms in the frame- 
work of Friedberg-Lee model to investigate the approach 
of strong interaction matter toward thermal equilibrium 
and the deconfinement phase transition during this pro- 
cess. The collision terms will be introduced through a 
relaxation time approximation in transport equations of 
quarks and sigmas, and the expansion of the system is 
simply described by the Bjorken scaling hydrodynamics 
[22]. We specially focus on the effect of collision terms 
and non-equilibrium on the deconfinement phase tran- 
sition, by investigating at what proper time the phase 
transition occurs, how long it lasts for a first-order tran- 
sition, and their dependence on the relaxation time. 

The outline of the paper is as follows. The full trans- 
port equations for quarks and sigmas in relaxation time 
approximation, and the simplified equations in quasipar- 
ticle limit and boost invariance approximation are pre- 
sented in Section 2. In section 3 we exhibit the numerical 
results and discussions. Finally a brief summary is given 
in the last section. 
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II. TRANSPORT EQUATIONS 



The Friedberg-Lee model is defined as [1-3] 



Cfl = i> {il^d^ - (mo + ga)) i> + -d'^ad^a - U{a), 



(1) 



where ip, 4> and a are quark, antiquark and scalar fields, 
respectively, TOq is the current quark mass and chosen to 
be in the following to simplify the calculations. There 
are five parameters in the Friedberg-Lee model, a with 
dimension h with dimension , dimensionless c, 
the coupling constant g between quark and scalar fields, 
and the bag constant B used to provide a quark confine- 
ment potential in mean field approximation in the vac- 
uum. Of the five parameters, two are adjusted to fit the 
proton size and the nucleon mass, and the others are left 
to survey, fit and predict physical data. Since the model 
is an effective; one;, the vahies of the parameters depend 
on the level of approximation employed [23]. Different 
parameter sets can be found in the book by Wilets [3] . 

The relativistically covariant quark Wigner operator 
Wq and sigma Wigner operator are the Fourier trans- 
form of the corresponding density matrices [15] 



d^ye'^^{x+^-)i>{x-^-) 



= I d'ye'pya{x+^)a{x-^) . 



(2) 



Calculating the first-order derivatives of $g and 
second-order derivatives of i>(7 with respect to x and y, 
and making use of the equations of motion 



(n^^M - {mo + ga))ip^ , 
d^'d^a + ^^+gU = 



(3) 



for the fields, one obtains evolution equations for the den- 
sity matrices. After Wigner transform one derives the 
kinetic equations (for details, see the similar work in [17] 
for QED) 



{K^K''-i:,)W,{x,p)=C, , 
{rK^-^q)Wq{x,p)=Cq 



(4) 



for the Wigner functions Wq{x,p) and Wa{x,p) which 
are, respectively, the ensemble average of the Wigner op- 
erators, 



Wq{x,p) = I d^ye'^{i>{x+^)i,ix-p , 
W,{x,p) = J d^ye^Py{a{x + |)a(:c - |)) . (5) 
The operators K^, and Eg are defined as 

Kt,=p^ + , 

E^, = (a + 6a + I ((a'a') + a^)) e'^'^^^^^S , 

Eg = (mo+gacos {^d^.d^j^ - ^ff^sm {^d^.d^^ , (6) 

where a and a' are, respectively, the mean field and quan- 
tum fluctuation of the scalar fleld, a = {a) and a = a+a'. 

The sigma and quark scalar densities {a' a') and {'tjj'tp) can 
be calculated through the Wigner functions 



< a'a' >= 



< i/jip >= Tr 



d^p 
(2^ 
d^p 



(27r) 



Wa{x,p) , 
Wq{x,p) . 



(7) 



In relaxation time approximation the collision terms C„ 



and Cq can be written as 



C„ = -ihp'^Uu, 



T th 



in ^ Wq-Wq*^ 



(8) 



where is the four-velocity of the hot medium formed 
by quarks and sigmas, and 6 is the relaxation time. It 
is necessary to note that the space-time derivative dx in 
the self-energies and Eg works only the mean field a 
and scalar density {a' a') on its left. 

Making ensemble average of the Klein-Gordon equa- 
tion in (3), one gets the equation of motion for the con- 
densate cr, 



(9) 



with the definition of the effective confinement potential 
Ueff{a, {a' a'), {'tpxjj)) at finite temperature 

dUeff dUia) {cr'a')„ , 

^^ = ^^ + ^Y^(& + ca)+5(V'^) (10) 

which depends on not only the mean field but also the 
sigma and quark scalar densities. The Eqs.(4) and (9) 
together determine the sigma and quark scalar distribu- 
tions and the condensate a self-consistently. 

From the first-order derivative of the effective potential 
f/e//, Eq.(lO), it can be defined as 
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U,ff{a,{a'a'),{i,i>)) = 

U{a) + U„{a, {a' a')) + Ug{a, {U)), 

Ua{a,{a'a')) = f dal-{a'a'){a,x) {b + ca) , 
Jo ^ 



(11) 



and Uq arise from the collective 



The two extra terms t/, 
motion of the quarks and sigmas. In the vacuum with- 
out collective motion the confinement potential Ueff = U 
has two minima, (Tper = corresponding to the pertur- 
bative vacuum and aphy 7^ to the physical vacuum, 
see the dashed line in Fig.l. Prom the definition of 
the bag constant, the energy density difference between 
the perturbative and physical vacua, B is determined by 
U{aphy) = 0, 



a 



B = — n-^ -I- — n-^ -I- — n-^ 



(12) 



The bag constant B depends strongly on the parameters 

a, 6, c and g. For various reasonable sets of parameters, 
however, B is around 20 MeV/fm'^ [3]. As discussed in 
[5] , it is the small bag constant B which leads to the low 
critical temperature of deconfinement phase transition in 
Friedberg-Lee model. 

An important aspect of the covariant kinetic theory is 
that the complex kinetic equation can be split up into a 
constraint and a transport equation [13,14,17 19], where 
the former is a quantum extension of the classical on-shell 
condition, and the latter is a covariant generalization of 
the Vlasov-Boltzmann equation. The complementarity 
of these two ingredients is essential for a physical under- 
standing of quantum kinetic theory. The quark Wigner 
function itself has no direct physical analogue since it is 
not a self-hermitian function. After the Lorentz decom- 
position [13,14,17-19], the kinetic equation for quark is 
changed into 16 transport equations and 16 constraint 
equations for the 16 Lorentz components of the quark 
Wigner function. In classical limit with ?i = 0, the con- 
straint equations for quark and sigma are reduced to 



[p" -ml)W„{x,p) = {) , 
{p'^ml)Wq(x,p) = Q 

with effective sigma mass and quark mass 



ml = a + ba +-{{&' a) +a'^ 



(13) 



= mo + ga 



(14) 



and the 16 spinor components are no longer fully inde- 
pendent. Only the quark number density fq and the 
spin density 50 are the fundamental elements and the 
other components can be expressed in terms of fq and 
go [17,18]. The classical transport equations for sigma 
and quark densities f„ and fq are reduced to the familiar 



Boltzmann equations in the rest frame of the heat bath 
where the four- velocity is = {1,0}, 



dtf<j,q + )/<T,g 



Z<i.v f 



fa,( 



fth 



with the particle energies Eg-^q 



m 



(15) 



The rela- 



tions between the scalar densities and number densities 
are 



m„ 



d^p 1 
(27r)3 
d^p 1 



(2^)3 Eq 



fa{x,p) , 
fqix,p) ■ 



(16) 



The equilibrium distribution functions in the classical 
transport equations (15) are the familiar Bose-Einstein 
and Fermi-Dirac distributions 



Tl 



(17) 



with the sigma and quark degenerates ga = ^ and gq = 
24. 

The classical transport equations can be greatly sim- 
plified by taking into account Bjorken's boost invariant 
picture often used to describe longitudinal expansion of 
relativistic heavy ion collisions where a central plateau of 
the final rapidity distribution exists [22]. Baym extended 
Bjorken's method of scalar hydrodynamics to the phase 
space and solved the Vlasov-Boltzmann equation with a 
constant particle mass [24]. 

Neglecting transverse expansion of the system and as- 
suming boost invariance along the longitudinal direction 
make the transport equations (15) for sigma and quark 
with effective masses much more simple, depending on 
the proper time r = 



only. 



th 



^rft7,q — 



Ja;q J a,' 



(18) 



with the solutions 

- ^ T 
/a,g(T,p) = e -^-0 "'■^ ^ fa,q{To,PT,Pz — ) 

To 



+ 



f 

J To 



0{r') 



f',%{T,pT,p.^) , (19) 



where the time ratios — and come from the longitu- 
dinal expansion of the system [24] . 

While the quark distribution in (19) is the same as 
that obtained in Baym scenario [24] , the quark mass here 
depends on the mean field a. This dependence couples 
the quark distribution fq with the meson distribution 
/ct and the mean field a. It is this coupling that leads 
to the phase transition from confinement to deconfine- 
ment through the change of the quark mass with the 
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mean field. In the longitudinal boost invariant picture, 
the equation (9) of motion for the scalar condensate a is 
simplified as a normal second-order derivative equation, 

(^52 + ldr^a+^ + ^ia'a'Kb + ca) + g{U) = • 

(20) 

The proper time dependence of the temperature, T{t), 
in the equilibrium distribution J*'' is determined by the 
energy conservation law in collisions, 

^^-^ [E^fa{T, (7,P)+ Egfg{T, (7,p)] , 

(21) 

where wc kept only the particle contribution to the en- 
ergy density, and neglected the mean field terms in e and 
e*'', since they are the same and do not affect the energy 
conservation law. 

In the relaxation time approximation, the key step is 
the computation of the relaxation time 9. In principle, 
depends on the type of particles, sigma or quark, and is 
a function of phase space coordinates. For all the discus- 
sions above we have neglected its type dependence and 
momentum dependence. Qualitatively, 9 has the order 
of the standard strong interaction scale, 6 ~ 1/m [24]. 
Considering the fact that collision terms are the driving 
force for the system to reach equilibrium, they affect the 
system strongly only in the beginning of the evolution 
and are damped with increasing time. With this feature 
in mind we model the relaxation time by a simple step 
function, 

where A is an inverse measure of the collision strength and 
Te indicates the collision duration scale. A has the order of 
1 fm and A ^ 00 means equilibrium without collision. In 
order to get a more realistic parametrization we replace 
in the following numerical calculations the function by 

a more gradual exponential function ^1 -|- e"^~T^^ with 
the collision duration variance /. When Z — > 0, it goes 
back to the step function. 

In order to solve numerically the transport equations 
for the distributions /^.g and the Klein-Gordon equation 
for the mean field a as functions of proper time r, one 
must know their initial values. Since the particles pro- 
duced in rclativistic heavy ion collisions at initial time 
To are essentially emitted from the colliding point z = 
at t = 0, the initial distributions in the central slice are 



peaked in the plane = [24] . Therefore, we can choose 
the initial distributions as 

/.,,(ro,p)=/*f,(To,p)<5(p,) . (23) 

With the known number distributions (19) one can get 
the scalar condensates {&'&') and (ipip) through the rela- 
tions (16) and in turn the effective quark potential Uejj. 

As we discussed above, there are two minima of the po- 
tential in the vacuum without collective motion, aphy is 
the global minimum but Gper is only the local minimum, 
see the dashed line in Fig.l. This means that the system 
is in physical vacuum. Since we focus in this paper on 
the phase transition from deconfinement to confinement, 
which can be realized in relativistic heavy ion collisions, 
we put initially the system in a deconfinement state with 
strong enough collective motion, and then study when 
the phase transition happens. To this end, the temper- 
ature To in the initial distributions (23) should be high 
enough to guarantee that the system is in perturbative 
vacuum initially, namely aper is the absolute minimum 
of the effective potential in the beginning. Due to the ex- 
pansion and the collision terms, the energy density and 
effective temperature of the system fall down, and the po- 
tential difference between the physical and perturbative 
vacua decreases during the evolution. At some critical 
time Tc the difference between the two vacua disappears 
and a first-order phase transition begins. 

We define Th as the terminating time of the first-order 
phase transition. During the process of the phase transi- 
tion, the deconfinement state and confinement state co- 
exist, the energy densities can be expressed as 

e(T, a) = x{T)e{T, aper) + (1 - x{T))eiT, (Jphy), 
e"'{Te,a) = x{T)e'\Te,aper) + (1 - x{T))e'\Te,cJphy), 

(24) 

where x(t) is the fraction of matter in the deconfinement 
phase at time Tc < t < th- Prom the energy conservation, 

e(T,a) =e*''(re,(7) , (25) 

one reads 
x{t) = 

e*''(Tc, aphy) - e(r, (Tphy) 

[e(T, aper) - e*'*(Tc, aper)] + [e*^{Te, aphy) - e(T, aphy)] 

(26) 

Clearly, at the beginning time Tc and the end time th of 
the first-order transition, one has 

x{Tc) = 1 , 

x{Th) = . (27) 
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III. NUMERICAL RESULTS AND DISCUSSIONS 

We choose the potential parameters [3,4,7] a = 
51.6/m-2 5 ^ -799.9/m-\c = 4000 and the coiiphng 
constant g — 14.8, which are satisfactory for fitting a part 
of the static properties of hadrons and lead to a small bag 
constant B = 2SMeV/ fm,^. The temperature in the ini- 
tial distributions (23) is taken as Tq = 160 MeV which is 
high enough to guarantee an initial deconfinement state. 
As explained above, the collision parameter A describes 
the intensity of collisions and therefore determines how 
long it needs to approach the equilibrium state, and Tg is 
the duration of the collisions which controls how long the 
interaction lasts. As discussed by Baym [24], A < 1/m 
is reasonable for strong interaction. According to the 
space-time scale of relativistic heavy-ion collisions, one 
can estimate Te ~ tqA^^^ fm. 
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FIG. 1. The effective potential Ueff as a function of a at 
initial time To, critical time Tc, and the limit r — > 00. The 
deconfinement phase transition begins at Tc = 5.6/m and ends 
at Th = 30 fm. 

Fig.l shows the effective potential Ueff as a function 
of a at different time for collision parameters A = 0.5 fm 
and Te = 4 fm. At initial time tq, Ueff has two min- 
ima, the local one at a = aphy and the global one at 
a = (jper = 0. The system stays initially in the perturba- 
tive vacuum. As time goes on, the system expands and 
its energy density drops down gradually. While the po- 
sitions of the two minima remain unchanged during the 
evolution, the potential difference between the two min- 
ima, the effective bag constant Beffir), becomes smaller 
and smaller. At a critical time Tc, Beff{Tc) = 0, the 
first-order phase transition begins and the confinement 
phase appears. After Tc, the temperature parameter in 
f^'^q{T, p) remains a constant Tc, the physical vacuum and 
the perturbative vacuum coexist, namely the deconfine- 
ment phase and confinement phase coexist until another 
critical time Th when the transition is totally completed 
and the system is purely in confinement state. For the 
collision parameters used, the numerical calculation gives 
Te = 5.6 fm, Th = 30.4 fm, and Tc = 94 MeV. After Th, 
the temperature parameter drops down again, and the 



global minimum is located at a = (Jphy When t — > 00, 
the effective potential Ueff approaches to U in the vac- 
uum, shown as the dashed line in Fig.l. Fig.2 shows the 
evolution of the temperature parameter before Tc- 
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FIG. 2. The temperatv/re parameter as a function of proper 
time before the phase transition. 

The two critical times, Tc and Th, strongly depend on 
the relaxation time parameters, A and Te- A system with 
strong interaction, namely small A, needs only a short 
time to approach equilibrium state, and the tempera- 
ture parameter drops down rapidly. Correspondingly, 
the phase transition occurs early and the duration of the 
first-order transition is short. The parameter Te controls 
also the values of Te and Th- A large Te means a long 
last of collisions, the system needs only a short time to 
get to the points where the phase transition begins and 
ends. If the duration of the collision Te is long enough, 
the system will approach equilibrium state eventually. If 
Te is too small, however, the system will not approach 
equilibrium state unless the collision strength is strong 
enough. Fig. 3 shows the effective potential as a function 
of (T in the limit of A ^ 00. In this case, the only driving 
force of the phase transition is the cooling of the mean 
field system, the dynamics is purely refiected in the effec- 
tive particle masses and niq. The transport equations 
(19) are reduced to 



/(7,g(T,p) = fa,q{'^0,PT,Pz — ) 

To 



(28) 



If we neglect the collective motion of a since it is heavy 

enough, and consider the fact that before the phase tran- 
sition quarks are massless in the deconfinement state, the 
energy conservation law in the limit of A — > 00 gives ap- 
proximately the time evolution of the temperature pa- 
rameter, T{t) ~ Tq {tq/t)^^^ - Only when the thermal 
distribution in (23) is taken as Boltzmann distribution, 
we have exactly. 



T(r) 



1/4 



(29) 



If we put the initial state in a 3-dimensional equilib- 
rium state and take the mass to be zero, the longitudinal 
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expansion of the system will force the state in nonequi- 
librium afterwards. The energy conservation in the limit 
of no collisions leads to 



T(7 



5(r) 




(30) 



Since the quark mass is taken as a constant, m = in 
this special case, the quark distribution is decoupled from 
the meson distribution and the mean field, the system is 
reduced to the case in Baym scenario. Therefore, the 
result (30) has a similar structure to eq. (19) of the work 
by Baym [24]. 

The comparison between (29) and (30) is shown in 
Fig. 4. Since the function ^(t) < 1, the temperature with 
initial equilibrium state is always lower than that with 
initial nonequilibrium state for any tim(^ r > rn. 
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FIG. 3. The effective potential Ueff as a function of a 
at initial time Tq, critical time r^, and the limit t ^ oo in 
the case without collisions. In this limit, the phase transition 
starts at Tc = Wfm and ends at Th = 841 fm. 



A = 0.5 fm and Tg = 4 fm and 89 MeV for A oo, the 
time structure of the phase transition is controlled by the 
collisions. The relaxation time dependence of the critical 
time Tc and the duration At = Tf, — scaled by their 
values in the limit of no collisions is shown in Figs. 5 and 
6. Both decrease monotonously with increasing collision 
strength and/or increasing collision duration. When the 
collisions are strong enough and the collision duration is 
long enough, the phase transition will start before the 
collisions cease. This is clearly reflected in Fig. 5 where 
Tc at A = 0.2 fm no more changes when Te > 4 fm. With 
increasing A the starting time Tc and the duration At 
will approach finally the limit values of no collisions, they 
converge at 10.12 fm and 831 fm, respectively. The speed 
of the convergence is governed by the collision duration. 
When the collision duration is short enough, the system 
is not sensitive to the collisions in the region of weak 
strength. This is the reason why Tc and At converge 
very fast with small Tg. 
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FIG. 5. The beginning time of first-order phase transition 
scaled by its limit value without collisions as a function of the 
collision strength for collision duration Te = 2,4,6 fm. 
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FIG. 4. The time evolution of the temperature parameter 
with initial nonequilibrium state (solid line) and with initial 
equilibrium state (dashed line) in the limit of no collisions. 

While the collisions do not change the critical value of 
the temperature parameter remarkably, it is 94 MeV for 
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FIG. 6. The duration of the first-order phase transition 
scaled by its limit value without collisions as a function of 
the strength for collision duration Te = 2, 4, 6 fm. 

Fig. 7 shows the temperature parameter as a function 
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of time for fixed collision duration Tg = 4 fm but differ- 
ent values of collision strength. Each line starts at initial 
time To = 1 fm and ends at the critical time Tc which 
depends on the collision strength. After Tc the temper- 
ature does not change until the second critical time Th- 
With increasing collision strength, namely decreasing A, 
the system cools down more and more fast, and the phase 
transition happens more and more early. 

In principle, only one of the two collision parame- 
ters A and Tg is free. A system with strong interac- 
tion needs a short time to approach equilibrium, and 
vice versa. While it is difficult to determine the rela- 
tion between the two parameters in relaxation time ap- 
proximation, we can obtain a constraint condition. For 
a system with fixed collision duration scale Tg, the col- 
lision strength should larger than the maximum value 
corresponding to the well-known Bjorken scaling solution 
T/Tq = (tq/t)^^^ [22]. Since A is the inverse measure of 
the collision strength, we have the constraint 

Amm < A . (31) 

For Te = 4/to, Xrnin IS approximately 0.3/m. When A is 
less than Xmin, the result will not be physical. 
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FIG. 7. The temperature parameter T as a function of 
proper time r at fixed Te ■ The solid curves from the bottom to 
the top correspond to \ = 0.5, 1 and 2 fm, and the two dashed 
lines correspond to the limit (29) and the Bjorken limit, re- 
spectively. 

IV. SUMMARY 

We investigated the deconfinement phase transition of 
an expanding quark system in Friedberg-Lee model, con- 
sidering the collision terms in relaxation time approxi- 
mation. We calculated the beginning time and the du- 
ration of the first-order phase transition for different col- 
lision strength. While the critical temperature of the 
phase transition is dominated by mean field, described 
by A ^ 00 in our treatment, the collision terms have sig- 
nificant influence on the beginning time and the duration 
of the transition. Strong collisions result in an early and 
short transition, and weak collisions make the transition 



begin late and last a long time. Although our results 
are derived in a particular model, we expect the quali- 
tative dependence of the phase transition times on the 
collisions to be of more general validity, because it is well 
known that collisions are the driving force of thermaliza- 
tion and control the speed of thermalization process of 
any system. 

Since the deconfinement phase is only an intermediate 
state in relativistic heavy ion collisions, there is no way 
to measure its properties directly. One needs to extract 
signatures of the phase transition and the new state from 
analysis of final state distributions of different particles, 
such as low moment dilepton enhancement, J/"^ suppres- 
sion and strangeness enhancement at SPS [25] and disap- 
pearance of back- to-back jets at RHIC [26]. Normally the 
theoretical study on these signatures is based on the as- 
sumption of equilibrium system. According to our discus- 
sion above, the thermalization process due to collisions 
among the particles make the phase transition happen 
early and shorten the duration of the first-order phase 
transition. As a consequence, the contribution from the 
deconfinement phase and the coexisting phase to the final 
state distributions will be changed, especially for those 
signatures extracted from the early evolution. 
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